Internal Friction and Vulnerability of Mixed Alkali Glasses 
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Based on a hopping model we show how the mixed alkali effect in glasses can be understood if only 
a small fraction cv of the available sites for the mobile ions is vacant. In particular, we reproduce 
the peculiar behavior of the internal friction and the steep fall ( "vulnerability" ) of the mobility of 
the majority ion upon small replacements by the minority ion. The single and mixed alkali internal 
friction peaks are caused by ion-vacancy and ion-ion exchange processes. If cy is small, they can 
become comparable in height even at small mixing ratios. The large vulnerability is explained by a 
trapping of vacancies induced by the minority ions. Reasonable choices of model parameters yield 
typical behaviors found in experiments. 

PACS numbers: 66.30.Dn,66.30.Hs 
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The mixed alkali effect (MAE) is a key problem for un- 
derstanding ion transport processes in glasses and refers 
to strong changes in transport properties upon mixing of 
two types of mobile ions (for reviews, see 0, 0)- Fun- 
damental for the MAE is the behavior of the tracer dif- 
fusion coefficients Da and Db of two types of ions A 
and B. With increasing number fraction x of B ions, 
i.e. with successive replacement of A by B ions, Da de- 
creases while Db increases. These changes are reflected 
in the activation energies Ea,b(x) ~ ~ k^T log D a, b(x), 
so that at low temperatures Da,b(%) vary by many orders 
of magnitude. As a consequence, the ionic conductivity 
a(x) ~ (1 — x)Da(x) + xDb(x) runs through a deep min- 
imum close to the intersection point of Da and Db- 
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ses of neutron and x-ray scattering data with the reverse 
Monte Carlo technique ^3 support the picture that a 
structural mismatch effect leads to distinct preferential 
diffusion pathways for each type of mobile ion. The con- 
nectivity and interference of these pathways determines 
the long-range ion mobilities. On this theoretical basis 
alone, however, main properties of the MAE are still not 
well understood. 

One property is the large "vulnerability" , that means 
the steep decrease of the diffusion coefficient of the major- 
ity ion, e. g. A, w ith beginning replacement by the minor- 
ity ion B [lfiUlfij. In many mixed alkali glasses, In Da (a;) 
is a convex function for small x, reflecting that Ea{x) is 
a concave function in these systems 0, 01 • Why is the 
increase of Ea largest at beginning replacement, where 
the influence of the minority ion should be weak? 

Another quantity poorly understood so far is the in- 
ternal friction (for a review, see 19]). Measurements on 
mixed alkali glasses show that the ionic motion leads to 
two mechanical loss peaks, a single alkali peak (SP) and 
a mixed alkali peak (MP). The peak frequency r" 1 of 
the SP has (nearly) the same activation energy E s as the 
diffusivity of the more mobile ion, while the lower peak 



frequency r" 1 of the MP exhibits an activation energy 
E m that is not simply related to that of the mobility 
of the less mobile ion. Particularly puzzling is that the 
height H m of the MP can be comparable to the height 
H s of the SP even at very small x, and that H m becomes 
larger the more similar the two types of ions are [2(1 ]. 

In this Letter we show that these peculiar behaviors 
can be explained, if the fraction cy of vacant sites for 
the mobile ions is small. Small values cy = 5 — 10% 
have been found recently in molecular dynamics simu- 
lations |2il |22| , and should be expected also on general 
grounds, since empty sites correspond to local structural 
configuration of high energy |23| . Here we argue that 
the experimental results for the internal friction provide 
independent support for this feature. 

To reason that conjecture, we note (see eq. (J3J) be- 
low) that the SP and MP can be attributed to processes, 
where the more mobile ions (e.g. A) exchange sites with 
vacancies (AV exchange) , and where the A ions exchange 
sites with B ions (AB exchange) . When denoting by ca,b 
the fractions of sites occupied by A,B ions, a mean-field 
argument then predicts that H s oc cacv and H m cx cacb 
should become comparable if cy ~ cb- Therefore, if it 
is possible to find H m ~ H s at small x, where cb ~ x is 
small, cy must be small also. Indeed, we will show in the 
following that this rough argument can be substantiated. 

First we are confronted with the general problem 
whether the mismatch concept is sustainable as mech- 
anism for the MAE if cy is small. To this end we inves- 
tigate a lattice model, where the rate wfj for a jump of 
an a ion (a = A, B) from site i to a nearest neighbor site 
j is determined by the site and barrier energies 
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The mismatch effect is taken into account via the mis- 
match energies ef > 0, u^ s > 0, and the structural 
variables (if, < < 1 that specify the a-character 
of site i. A site i with large fif has a local environment 
favorably accommodated to an a ion and hence a low 
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energy ef. Large values of /if imply a small value of 
/if and vice versa. This is accounted for by taking op- 
posing mean values /if — 1 — /if in the distributions of 
the /if (see below). Below a threshold value /if, sites 
loose their identity and we require the energy to saturate 
(ef = -e mi? /if for /if < /ifj. 

The barriers itfj- in eq. Q contain a bare structural bar- 
rier Uq (characterizing the activation energy in the pure 
systems) and a mismatch barrier that becomes higher 
with increasing "foreign" (3(^ a)-characters of the initial 
and target sites involved in a jump of an a ion. This is 
necessary in order to obtain low mobilities for the minor- 
ity ion in the dilute regime |2*3 |. Within an Anderson- 
Stuart like picture 25] it can be physically understood 
by the fact that larger minority ions have to bring up 
an additional elastic strain energy when open up door- 
ways, and that smaller minority ions have to surmount a 
higher saddle point energy associated with the Coulomb 
potential of the counter ions. 

The disorder in the glass is reflected in the distribu- 
tions of the /if . Due to the local accommodation of 
the network structure to the ions in course of the freez- 
ing process, the numbers of sites with large /if and /tf 
scale with (1 — x) and x, respectively. In addition we 
take into account short-range correlation between the /if. 
A convenient (technical) way to generate corresponding 
distributions (the detailed form is not important here) 
is as follows: We draw /if from truncated Gaussians 
(p(/if) oc exp[-(/if - /2f) 2 /2A^] for fxf e [0,1]) with 
/if = Vj/( z + 1) = 1 — A*f j where the sum runs over 
site i and its nearest neighbors, and r\j are random binary 
numbers equal to one with probability (1 — x). 

Kinetic Monte Carlo simulations are performed 
in a simple cubic lattice with periodic boundary 
conditions, using Metropolis hopping rates wfj = 
^ Q exp(-uf J /fc B T)min(l,exp[-(e° ! - ef)/k B T}) {v a are 
attempt frequencies). We choose cy = 0.05 and con- 
sider, for computational convenience, a symmetric set of 
parameters v K = v^ = v, ef is = ef lis = e mis , ... To deter- 
mine the tracer diffusion coefficients and internal friction 
as functions of x and reduced temperature 9 = k B T /e m i s , 
we fix the remaining parameters to /i c = 2/7, A = 1/14, 
and M m i s = 0.7e m ; s . The parameter e m j s characterizes the 
difference between the two types of mobile ions (e.g. with 
respect to size). 

Figure shows the normalized diffusion coefficient 
Da (x) /Da (0) as a function of x for 5 different 9 [due 
to the symmetric parameters we do not plot D B (x) — 
Da(1 — x)). The curves exhibit the typical behavior found 
in experiment: Da(x) decreases by many orders of mag- 
nitude when A ions are replaced by B ions. This effect 
becomes stronger with lower temperature or stronger dif- 
ference e m i s between the two types of mobile ions, i.e. 
with larger 9 = k-QT/e m { s . We thus conclude that the 
mismatch concept is sustainable for small cy. 

In particular, \og[DA(x)/ Da{0)] = —EA{x)/k B T, ex- 




FIG. 1: (a) Normalized diffusion coefficient Da{x) / Da{0) of 
A ions for 5 different reduced temperatures 6 = k^T/t-mis. 
The inset shows Da{x) for small x at 6 = 1/10 (open sym- 
bols) in comparison with the results in the absence of a cor- 
relation induced trapping effect (full symbols), (b) Measured 
data (symbols) for (1 — x)Na20-a;Rb2-4B203 glasses at 652K 
(redrawn from |18|1 and fit by the model (lines) with = 1/9 
and asymmetric mismatch barriers «^ is = 2uf is = e m i s . 

hibits a convex curvature for small x, and Ea(x) a con- 
cave curvature (see Fig.0Ji; note that one must add Uq to 
obtain the total activation energy E A ot (x) = Ea(x)+uq). 
The large vulnerability is caused by a trapping effect: For 
small x, there exist regions in the glass where the network 
structure tended to accommodate to the "foreign" B ions. 
Because of the spatial correlations, the regions consist of 
several sites with smaller A character /if than on aver- 
age. These sites are not only favorably occupied by B 
ions, but also favorably occupied by vacancies, because 
this allows more of the majority A ions to occupy well ac- 
commodated sites with large /if close to one. Vacancies 
hence become trapped in regions more favourably acco- 
modated to the the foreign B ions and they have to bring 
up an additional activation energy to promote the diffu- 
sion of A ions. The corresponding reduction of Da(x) 
is strongest for x — ► 0, since with increasing x the num- 
ber of sites with values /if close to one decreases. The 
importance of the trapping effect |26| is demonstrated in 
the inset of Fig.^i, where Da(x) is shown in comparison 
with the curve obtained in the absence of spatial corre- 
lations (/if = ryf instead of /tf = ^jVj/i 2 + 1)) but 
otherwise the same parameters). 

Fig. resembles well the experimental behavior. By 
choosing asymmetric parameters one can reproduce mea- 
sured curves as shown in Fig. \I]p with reasonable param- 
eter values (if we require 9 = 1/9 in Fig.^Jj to represent 
the experimental temperature 652K, we would obtain a 
mismatch energy e m ; s = 0.5eV). A detailed fitting to ex- 
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FIG. 2: (a) Internal Friction Q and spectral components 
(cf. eq. |5J as a function of u for x = 0.3 and 9 = 1/10. (b) 
Change of the internal friction Q _1 with x for 9 — 1/10. 

perimental results, however, is not our aim here (and has 
limitations due the simplicity of the model). 

Rather, we next explore how the behavior in the inter- 
nal friction Q^ 1 (u>, T) can be understood. To this end, we 
consider a small modulation Aef(t) = ao K[x"e _Mt ] of 
the site energies in response to an oscillatory shear stress 
with amplitude ao and frequency u> |27|. The complex 
coupling constants xf fluctuate due to the disorder in the 
glass. Without making specific assumptions regarding 
their correlation properties, we employ a random phase 
ansatz and regard them to be uncorrelated for different 
sites i and different types a, [xf'Xj* ]av °c and 
uncorrelated with the site energies of the unperturbed 
system [x? e i*]av = ([. . .] av denotes a disorder average). 
Introducing the occupation numbers nf(t), i.e. nf(t) = 1 
if site i is occupied by species a (a = A,B,V here) at 
time t and zero else, linear response theory then expresses 
the internal friction in terms of the Fourier cosine trans- 
forms of the exchange correlation functions (nf(t)rij(0)) 
((. . .) denotes a thermal average). The result is 

Q~ l {u>,T) = c aV S a y + 2c AB S AB \ (2) 

B a— A,B 
/>oc 

S aP (u>,T) = / dtC a0 (t)cos(ojt), (3) 
Jo 

where n is constant 0, c a p — [(nf )(nf )] av , and C a/ 3(t) = 
1 — i( n f W n f }]av/ca/3 are the normalized exchange corre- 
lation functions \C a fj(Q) = 1 and C a p{t — > oo) = 0]. 

Figure^ shows Q~ x as a function of u> for x = 0.3 
and 9 = 1/10, together with the spectral components 
Q~p oc (u>/ 9) c a pS ' a p in eq. J3J). Two arise from exchange 
processes of the ions with the vacancies (AV and BV 
exchange) and one arises from exchange processes of the 
two types of ions ( AB exchange) . The exchange of vacan- 
cies with the majority ions (AV) leads to the SP and the 
AB exchange to the MP, while the BV exchange yields 
a broad spectral contribution that is masked by the SP 
and MP, and cannot be resolved in The overall be- 

havior of Q^ 1 with varying x for 9 = 1/10 is shown in 
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FIG. 3: (a) H m in dependence of # _1 = emis/hsT at x = 0.5. 
The inset shows the decrease of cab with increasing 9~ 1 oc 
£mis (from analytical calculation), (b) _ff m as a function of 
ion radii ratio of various mixed alkali pairs for x — 0.5 and 
lo/2-it = 0.4Hz fixed. 

Fig. 03. Due to the small c v , the MP can be well identi- 
fied already for x — 0.01, in agreement with the argument 
given above. With increasing x the MP becomes higher 
and moves towards larger frequencies, while the SP be- 
comes lower and moves towards smaller frequencies. At 
x = 0.5 the SP cannot be resolved any longer. 

Remarkably, also the puzzling rise of the height i/ m 
with increasing similarity of the two types of mobile ions 
(i.e. decreasing e mis ) is reproduced, see Fig. H m at 
x = 0.5 becomes higher with decreasing e m i s (we regard T 
to be fixed here). The reason for this behavior lies in the 
weighting of the AB peak by the disorder averaged prod- 
uct cab = [{n 
(nt), (nf) in eq. © 
and B ions share the same sites with lower probability 
(see inset of Fig. an d accordingly H m decreases. The 
behavior is reminiscent to the variation of H m with the 
fraction of ion radii observed in experiments at x = 0.5 
for a fixed frequency, see Fig. 13b |2j|. 

When denoting by E a p the activation energies of the 
Q~a peak frequencies , we can define by E s — 
min(i?AV! -Ebv) an d E m = Eab the activation energies 
of the SP and MP. The activation energies E a p obtained 
from Arrhenius plots (for the same 9 values as in Fig. ^) 
are shown in Fig. 0Ji together with the diffusion activa- 
tion energies E < = mm(E A , E-q) and E > = max(i?A) £-b) 
of the more and less mobile ion. We find that E s 
(= min(SAV) Ebv)) follows closely E <7 in agreement 
with measurements 0. By contrast, E m (= E AB ) dif- 
fers from E > . The deviations found in Fig.0Ji are similar 
to those found in experiments, see Fig.0Ja |3J|. 

To show Q^ 1 as function of temperature for fixed fre- 
quency, as typically obtained in experiments, we use 
time-temperature scaling to transform Q _1 from fre- 
quency to temperature space. Taking the 9 values from 
Fig. ffli, we first checked, that a scaling ujS a p(uj,T) — 
Faf}{wr a p(T)) is approximately valid. By applying this 
scaling we then replot Q _1 from Fig. [2b as function of T 
in Fig. 0t for small x and parameter values e m i s = leV, 
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FIG. 4: (a) Activation energies of the tracer diffusion coef- 
ficients shown in Fig. and of the spectral components of 
the internal friction, (b) En&, E^b, and E m in (1 — x)Na20- 
xRb20-3Si02 glasses (redrawn from fl^l '). (c) Q~ l calcu- 
lated from Fig. |2] after applying time-temperature scaling for 
e m i s = 2lto = leV and lo/v = 10 12 (see text), (d) Q^ 1 in 
(l-a;)Na20-a;Rb20-3Si02 glasses for lu/2ty = 0.4Hz (redrawn 
from Q). 

uq = 0.5eV, and uj/v — 10~ 12 . These values yield activa- 
tion energies E t ^ t {x) = Ea + u ~ 0.6 . . . 1.2eV compara- 
ble to those found in experiments, cf. Fig.^J). The choice 
uj/v = 10~ 12 corresponds to ui ~ 1Hz, when v ~ 10 12 Hz 
is a typical attempt frequency. Using this rough param- 
eter estimate, peak temperatures in Fig. 0J; are obtained 
that reflect the experimental scenario in Fig. Hi [note 
that the glass transition and the "non-bridgening oxy- 
gen" (nbO) peak 0] are not taken into account in the 
present approach]. 

In summary, we have shown how the mixed alkali ef- 
fect in glasses, including the vulnerability and the pecu- 
liar features of the internal friction, can be understood 
based on the mismatch effect in the presence of a small 
vacancy concentration. Due to the small cy, H m can 
become comparable to H s even at small concentrations 
of the minority ion, and H m increases, when it becomes 
easier for both types of ions to share the same sites. The 
large vulnerability can be connected to a trapping of va- 
cancies induced by the minority ions. Reasonable choices 
of parameters allowed us to faithfully reproduce typical 
behaviors found in experiments. 
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